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Abstract. 

' r- '. We investigate the time evolution of three recently introduced distance-type 

[ measures of non-Gaussianity for one-mode states of a radiation field in contact with a 

thermal reservoir. They are based on the Hilbert-Schmidt metric, relative entropy 
and Bures metric. We first write the damped density matrix for an input Fock- 
diagonal state using the Weyl expansion of the density operator. Specifically, for 
an M-photon added thermal state, we obtain a compact analytic formula of the 
time-dependent density matrix that is used to evaluate and compare the three non- 
Gaussianity measures. We find a good consistency of the examined measures on the 
sets of damped states. Finally, we show that gaussification of the attenuated field 
mode is accompanied by a non-monotonic evolution of the von Neumann entropy of 
its state conditioned by the initial value of the mean photon number. 



PACS numbers: 42.50.Dv; 42.50.Ex; 03.65.Yz 
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1. Introduction 

In quantum optics a non-classical state is denned as having a negative or a highly 
singular P representation (more singular than Dirac's 5). Otherwise we term it as a 
classical state. Cahill p] and later on Hillery j2j proved that the only pure states that 
are classical are the coherent ones: all other classical states are mixtures. On the other 
hand, according to Hudson's theorem [3], all pure states except for the Gaussian ones 
have negative Wigner functions. We can thus remark that all pure states with negative 
Wigner function are both non-classical and non-Gaussian. In such cases, negativity of 
the Wigner function could be interpreted as an indicator of non-classicality as prominent 
as the negativity of the P representation. Moreover, signatures of non-classicality 
were identified through the negativity of the Wigner function for some mixed non- 
Gaussian states as well [I]. Interest in the non-Gaussian states has recently emerged 
in quantum information processing. It was realized that non-Gaussian resources and 
operations could be more performant in some quantum protocols such as teleportation 
[5j [6j [7] and cloning [8] or even indispensable in entanglement distillation |9J. To 
quantify the non-Gaussianity as a resource in such cases, some distance-type measures 
of this property were proposed PU EH] • However, in order to avoid the complications 
encountered in extremization procedures of the previously defined distance-type degrees 
of entanglement p2], another pattern was followed in the case of non-Gaussianity degrees 
[T0| [TT] . Let us denote by pc the Gaussian state having the same average displacement 
and covariance matrix as the given state p. This state was reasonably used as a reference 
one in defining a degree of non-Gaussianity as the distance between the given state and 
the Gaussian state pc- In Refs. [TO f [TT ] IT3] , Genoni et al. employed the Hilbert-Schmidt 
metric to define a corresponding non-Gaussianity measure, 



and the relative entropy to introduce the entropic measure of non-Gaussianity: 



On the one hand, Eq. (11. ip gives an easily computable expression. On the other hand, 
despite its not being a true distance, the relative entropy is acceptable and used as a 
measure of distinguishability between two quantum states. It is remarkable that the 
relative entropy (11.21) reduces to a difference of von Neumann entropies: 



where S(p) := — Tr(plnp) is the von Neumann entropy of the state p. General 
properties of the non-Gaussianity measures (II. ip and (11.31) were discussed in detail 
in Refs.JJUlITTjITS]. Another approach to quantifying non-Gaussianity is based on the Q 
function and leads to a measure expressed by the difference between the Wehrl entropies 
of the Gaussian state pc and the given non-Gaussian state p [14J. It is also worth 
noting that the non-Gaussianity measures [10[ [11] were used to check on the relation 
between non-Gaussianity and the possibility of extending Hudson's theorem to mixed 
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states [H)J [IE]. Interestingly, non-Gaussianity in terms of relative entropy [TT] was 
experimentally measured for single-photon added coherent states [17]. Very recently, 
an experimental study on the same measure was reported for phase-averaged coherent 
states [T8] . 

The aim of this paper is threefold. First, we reinforce the Bures distance to the 
Gaussian state pc as a measure of non-Gaussianity of an arbitrary one-mode state p. We 
have recently introduced it on general grounds in Ref . [19] by using a well-known metric 
related to the fidelity between two quantum states [20] . Our motivation to introduce this 
measure was given by the valuable distinguishability properties of the fidelity described 
in Ref. [211 122]. Application of this measure for M-photon added thermal states showed 
us its good agreement with previously defined distance- type measures [101 [TT1 [13] . Note 
that fidelity-based metrics have proven to be fruitful in quantum optics and quantum 
information as measures of nonclassicality [23] and entanglement [21 1 [25 | 126]. Second, 
we use the three above-mentioned distance-type measures to investigate the loss of non- 
Gaussianity for M-photon added thermal states of a field coupled to a heat bath. To this 
end, we give a new derivation of the Fock-basis solution of the quantum optical master 
equation by employing very conveniently the Weyl expansion of the density operator. 
Third, we examine the possible non-monotonic evolution of the von Neumann entropy 
during decoherence. 

The paper is organized as follows. In Sec. 2 we recapitulate the principal features of 
the Bures degree of non-Gaussianity arising from some general properties of the fidelity. 
We here concentrate on an easily workable case, namely, that of mixed states having 
diagonal density operator in the Fock basis. A comparison between various degrees of 
non-Gaussianity is further intended. Therefore, we consider in Sec. 3 a set of states 
generated as solutions of the quantum optical master equation [2?]. The density matrix 
of the damped field at the moment t is derived for an arbitrary input state which is 
diagonal in the Fock basis. Section 4 is dedicated to a specific damped state which 
is important in applications: an M-photon added thermal state evolving under the 
influence of a thermal reservoir. With the explicit input density matrix at the moment 
t = written in Ref. [19] and exploiting some results of Sec. 3, we here find the damped 
density matrix in a closed analytic form. Further, the time evolutions of the Hilbert- 
Schmidt, entropic, and Bures measures of non-Gaussianity are analyzed and compared. 
A special attention is paid to the entropy production S(p) as another tool to examine 
gaussification of single-mode states due to the field-reservoir interaction. Section 5 
presents a discussion of our results focused on the consistency of the analyzed measures 
of non-Gaussianity. 

2. Bures measure of non-Gaussianity 

Let us consider two arbitrary states, p and a, of a given quantum system. According 
to Uhlmann, when both states are mixed, a good measure of the closeness between 
their properties is the maximal quantum-mechanical transition probability between 
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their purifications in an enlarged Hilbert space [20J [21]. This is an extended notion 
of transition probability between quantum states which is now called fidelity [21] and 
has the intrinsic expression [20] 




In Ref.[19] we have defined a fidelity-based degree of non-Gaussianity 

MP] := \ Pa) = 1 - VF(pJg). (2.2) 

As seen in Eq. (12.21) the fidelity is tightly related to the Bures metric ds introduced in 
Ref. [28J on mathematical grounds. Several general properties of fidelity are listed below 
together with the features they transfer to our definition (12. 2p . 

PI) < F(p, cr) < 1, and J-"(p, a) = 1 if and only if p = a. This property implies: 

5f\P) =0, iff p is Gaussian, (2.3) 
< 5 F [p] < 1, otherwise. (2.4) 

Properties ( 12. 3 p and (12. 4p justify the interpretation of 5p[p\ as a degree of non- 
Gaussianity. 

P2) F(a, p) = F(p, a), (symmetry). 

P3) J-"(p, a) > Tr(pcr); when at least one of the states is pure, then the inequality 
is saturated: ^(a, p) = Tr(pa). Accordingly, for any pure state |^ r )(^ r | Eq. ( 12. 2p 
reads: 

(J,.[|*)(*|] = 1- VWd*>- (2-5) 

Unlike 5 re which has the same value for all pure states possessing the same 
covariance matrix, 5f is state-dependent. 
P4) Fijj pW ,UaW) = ^(p, cr), (invariance under unitary transformations). As 
shown in Refs. [TO] fTTl [13], when U are the unitary operators of the metaplectic 
representation on the Hilbert space of states we have that 

p> = u0 (p') G = ijfcifl 

and therefore 

6 F [Upfr] = 6 F \p]. (2.6) 

It follows that, in the one-mode case, Sf[p] does not depend on one-mode squeezing 
and displacement operations. 
P5) J-(<&(p), > J-{p,cY) (monotonicity under any trace-preserving, convex- linear, 

and completely positive map $). This means that fidelity does not decrease 
strictly under any trace-preserving quantum operation, including non-selective 
measurement. Therefore, by virtue of the definition ( I2.2p . the Bures degree of 
non-Gaussianity does not increase under any such quantum operation. 
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P6) J-{f>i ® p2,o~i g) cr 2 ) = ^(pi, <3"i)J-'(p2, <t 2 ), (multiplicativity). Let us consider a 
bipartite product state p := pi ® p 2 . If p 2 is a Gaussian state we get 

F{pi ® P2, (pi)c ® P2) = (pi) G ) 

and therefore 

<Wi®P2] = 5 F [px\. (2.7) 

P7) For commuting density operators, [p, a] = 0, Eq. ( 12. ip simplifies to 

^(p,a) = [Tr(v/pV^)] 2 . (2.8) 

Properties (12. 3p . (12. 6p . and (12. 7p of 5p[p\ are shared by the non-Gaussianity measures 
( II. ip and ( II .3p as well p2]. In general, the Hilbert-Schmidt measure ( ll.ip is easier 
to compute than 6~f[p] ot 5re[p] and one could wonder why these last two were in 
fact considered? The answer is given by the monotonicity property under any trace- 
preserving quantum operation shared by the relative entropy and fidelity, unlike the 
Hilbert-Schmidt distance [22] . One of our goals in this work is to notice the role the 
monotonicity plays in evolution of non-Gaussianity under damping (see Sec 3). 

In the following we concentrate on an interesting computable case. The Gaussian 

reference state of a Fock-diagonal state, 

00 00 

p = J2pi\l)(l\ with X> = 1, (2.9) 
1=0 1=0 

is a thermal state with the same mean occupancy (N) = ^ z Ipi, 

PG=?s l \l)(l\ with si = — ^—l . (2.10) 

((Ao + iy +i 

The corresponding Hilbert-Schmidt and entropic degrees of non-Gaussianity were 
written in Refs. [101 [11] as: 



<Wp1 = 2 



1 + E^(^TI- 2 ^ 



(2.11) 



where Tr(p^) = (2(N) + l)" 1 is the purity of the state pc, and 

<Mp] = J> ln P l + + X ) ln ^> + X ) - ft) In(iV). (2.12) 

1 

In the last line we have used the von Neumann entropy of a Gaussian state. In this 
special case we notice the commutation relation [p, p G ] = 6 which allows the use of Eq. 
(£3D to get 



DC 



8 F [p] = 1-J2VpT^i. (2.13) 
1=0 

Note that various excitations on a thermal state of the type p ~ (a^) k a 1 p t h (&) a k 
have the density operator of the type (12.91) . Here a and are the amplitude operators 
of the field mode. 
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3. Master equation of one-mode field damping 

Loss of non-classical properties of the single-mode field states (Gaussian and non- 
Gaussian as well) during the interaction with a dissipative environment has been 
intensively studied in the last decades [29j [30j EU E2J EHJ [34]. Quite recently, more 
general properties of damped non-Gaussian states such as the evolution of mixing 
measured by linear entropy production [351 ESI EI] , or the evolution of various measures 
of non-classicality [38J were investigated by employing the quantum optical master 
equation in the interaction picture: 

= — (2apa^ — a)ap — pa) a) + jnn(a)pa + apa) — a)ap — pad)). (3.1) 

In Eq. (13. ip . p is the reduced density operator of the field, 7 is the coupling constant 
between field and bath, and Ur stands for the mean occupancy of the reservoir. This 
popular master equation is of the Lindblad type and thus preserves both the positivity 
and normalization of the density operator. It has the clear physical significance of 
describing decoherence of a field mode coupled to a heat bath. For a comprehensive list 
of references regarding solutions arising from Eq. (13. ip for the expectation values of field 
operators, characteristic functions and quasiprobability densities, we refer the reader to 
some recent work of Dodonov [39l HQ]. We here intend to write the evolving damped- 
mode state generated by this master equation from an input non-Gaussian one. Then 
we will compare the time evolution of the above-discussed degrees of non-Gaussianity. 
To go on this programme, we need a workable solution of Eq. (13. ip in the Fock basis. 
In what follows our main tool is the Weyl expansion of the density operator: 

p= l -J d 2 A x(A)D(-A). (3.2) 

In Eq. ( 13. 2ft . x(A) is the characteristic function (CF) of the state p, defined as the 
expectation value of the displacement operator D(X) = exp(Ad^ — A*d), 

X (X) := Tr[pD(X)]. (3.3) 

The CF x(X,t) of the damped field state is found to be determined by its initial form 

x(A,o) BDEH1EZ]: 

X (A,t) = x(Ae-i',0)exp - U R + (1 - e^)|A| 2 . (3.4) 

The mean photon number in the damped mode is 

(N(t)) = (N(0))e-^ + n T (3.5) 

where ny := tir(1 — e -7 *) is the thermal mean occupancy in the field mode at time t 
and (N(0)) is the initial mean photon number. 

We expect that an input Fock-diagonal density operator remains diagonal under 
the damping master equation (13.11) which is phase-insensitive. Indeed, the CF of an 
input Fock-diagonal state p is 

/ |\|2\ 00 

X(A,0) =exp J>(0) M|A| 2 ), (3-6) 
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where pu(0) is the photon- number distribution and Li is a Laguerre polynomial of degree 
I. Equation (13. 4p gives further 



x(A,t) = exp 



n T + - ) |A| 



X>(0) ^(|A|V 



(3.7) 



The density matrix Pjfc(t) of the damped field state is obtained from the Weyl expansion 
( 13. 2 p after here inserting the CF f !3.Tj) . Use is also made by the matrix elements of 
the displacement operator in Fock basis. After a simple calculation using the polar 
coordinates and an obvious change of variable in the integral (I3.2p we are left with the 
following expansion 



Pjk(t) — 8jk 



°° POO 

^Ph(O) / dx exp [-(fir + l)x] ^(xe -7 *) L 3 



(3.8) 



The integral in the above equation can be routinely performed yielding the formula 



Pjk(t) — 8jk- 



(n T y 



1=0 



[n R + l)(l-e 
n T + 1 



x 2 F X -j, — 



-it 



, , • (3.9) 

where 2 -F\ is a Gauss hypergeometric function, Eq. f ]A.ip . Equation ( 13. 9 p is the general 
solution for the damped density matrix as function of the input one in the Fock-diagonal 
case. Let us take the limit t — > oo in Eq. ( 13. 9p . The result is a thermal state with mean 
occupancy n R . We thus deal with a gaussification process, namely, evolution under 
this master equation eventually destroys the non-Gaussianity and also non-classicality 
properties of any input state. A special case of Eq. (I3.9P arises for the thermal contact 
with a zero-temperature bath [n R = 0), i.e., for the field coupling to the vacuum. A 
more familiar result is then obtained: 



Pjk(t) 




Pll (0)v j (l-V) 



1-3 



(3.10) 



where r] := e -7 *. Equation ( I3.10p describes both dissipation by contact with a zero- 
temperature reservoir and photon counting for which 1 — 77 is called quantum efficiency 
of the detector 



4. A case-study: photon-added thermal states 

According to an interesting remark of Lee [12], adding photons to a classical state implies 
removal of the vacuum state from its Fock expansion. A nonclassical output is obtained 
in every case. Addition of photons to a classical Gaussian state, such as a coherent or 
a thermal one, generates a non-classical output which is no longer Gaussian. We can 
say that photon-added states are both non-classical and non-Gaussian. An M-photon 
added thermal state [131 HI] is an interesting example of a Fock-diagonal state whose 
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non-classicality was recently investigated experimentally [45, 46J. The density operator 
of the state obtained by adding M photons to a thermal one p t h is 

n( M ) - - (a r ) M on, a M (A I) 

The spectral decomposition of a given thermal state p t h is of the type (I2.1(jp . with a mean 
photon number denoted n. The non-classicality of the state (14. ip was first discussed 
by Agarwal and Tara [43J in terms of its non-positive P representation and Mandel's 
Q-factor. 

4-1. Non-Gaussianity of photon- added thermal states 

The non-Gaussianity of the state (14. ID was recently evaluated in Ref.[T4"] by employing 
the Wehrl-entropy measure and found to coincide with the non-Gaussianity of the 
number state \M)(M\, being thus independent of the thermal mean occupancy n. We 
employ here Eq. (14. ip to write down the density matrix: 

n l ~ M 

■ Pl$lm- (4-2) 




(n + 1)' +1 

The photon-number distribution of the associated thermal state, 

is obtained by using the shared mean occupancy (N) = n(M + 1) + M. The non- 
Gaussianity degrees (12. lip . ( 12.121) and ( 12.131) for the state (14. 1 1) were carefully analyzed in 
our work [19] . We have there adopted as a reasonable criterion for their appropriateness 
a monotonic behaviour with respect to the average photon number (N) of the state. In 
the case of Hilbert-Schmidt degree (12.111) we were able to give an analytic result while for 
fidelity-based degree and relative-entropy measure we performed numerical evaluations. 
We found that they depend on the thermal mean occupancy n unlike the Wehrl-entropy 
measure [14] . Our evaluations have also shown a consistent relation between these three 
non-Gaussianity measures. Their behaviour under damping is analyzed in the following 
by employing the results reported in Ref . |19j . 



4-2. Gaussification of photon- added thermal states 

Insertion of the input photon- number distribution, Eq. ( 14. 2[) . into Eq. ( 13.91) gives the 
Fock-diagonal density matrix of a damped photon-added thermal state. Technically, we 
have to evaluate the sum of a power series involving Gauss hypergeometric functions 
which is discussed in Appendix A. We apply Eq. (1A.5I) and get 



jj 



(w fl + l)(l-e 
n e -7 * + fiT + 1 



M 



(ne 7 * + n T ) j 
(n e-T* + n T + 



, (fi _|_ \) e ~7* \ 



(n R + 1)(1 -e" 7 *)(ne- 7i + n T ). 




yt yt 

Figure 1. Time evolution of the fidelity-based non-Gaussianity, Eq. (|2.f 3j) (left plot) 
and the entropic measure, Eq. (12.12)) (right plot), for input Af -photon added thermal 
states with M = 1, 3, 5, 10 (from bottom to top) and n = 1 in contact with a heat bath 
having hr = 0.1. 

The Gaussian state associated to this density matrix is a thermal one with the mean 
occupancy given by the average photon number of the damped field state. According 
to Eq. (13. 5p . we have 

(N(t)) = [n(M+ 1) + M]e" 7 * + n T . (4.5) 

We have used the time-dependent density matrix (14. 4 ft and the mean occupancy (14. 5 j) 
to evaluate numerically two of the distance-type measures of non-Gaussianity we are 
interested in: 5p[ffi M >\ and SreIp^} via Eqs. ( 12. 13ft and ( 12. 121) . respectively. Our results 
are displayed in Figs.l and 2 where the time evolution of these non-Gaussianity measures 
are presented for weak damping (ur = 0.1) and a large thermal mean occupancy 
(ur = 5), respectively. 5f\P^] and 8re[p^ m >\ are monotonically decreasing in time 
in accordance to the corresponding property discussed in Sec. 2. The two measures 
considered above are also very much in agreement when the reservoir is much noisier 
as in Fig.2. They show a good consistency by inducing the same ordering of non- 
Gaussianity when considering specific sets of damped states continuously generated by 
the field-reservoir interaction. The last degree of non-Gaussianity we are interested in 




yt yt 

Figure 2. As in Fig.l, but for input M-photon added thermal states with M — 
1, 3, 5, 10 (from bottom to top) and n = 0.5 in contact with a heat bath having fin. = 5. 
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is evaluated using Eq. (12. lip . Fortunately, as in the input-state case [JJJ], we are able 
to give nice analytic expressions for both the purity Tr { [p^ M ^(i)] 2 } of the damped state 
and the scalar product Tr[p^ M \t)pQ (t)]. We first write 



Tr 



W)) 



oo 

E 

j=0 



( (N(t)) 



(M) 
Pji 



(t), 



(4.6) 



where the expectation value (N(t)) is given in Eq. ( 14. 5ft . By inserting Eq. ( 14.41) for 
the damped density matrix into Eq. (14.61) we eventually get a sum of hypergeometric 

^3l) . After a little algebra we obtain the formula 

((N(t)) + nr + l-e-^) M 



functions of the type 



Tr 



P 



((N(t)) +n T + ne-^ + l) 



M+l 



Now, the purity of the damped state, 



oo „ 

^{[p (m) (0] 2 } = E[^w] 



(4.7) 



(4i 



3=0 



is proportional to the sum of a power series whose coefficients are squared Gauss 
hypergeometric polynomials. We apply the summation formula (1A.7j) and find after 
some minor rearrangements: 

[e-T*(ii + 1)] 2M 



[2n T + 2ne-7* + 1] 2M+1 

[2n T + ne'^ + 1 



-7*1 



x 2 F l ( -M,-M;l; 



(4.9) 



[e-7*(n + l)] 2 

Note that we can express the r.h.s of Eq. (14. 9p in terms of a Legendre polynomial via 
Eq. (IA.2I) . The equilibrium values of both the scalar product ( 14 ,7p and the purity (14.91) 
are equal to the purity of the thermal state eventually imposed by the reservoir, namely, 
(2fiR + l) -1 . The interesting case of a number state \M)(M\ is given by setting n = 
in the above formulae. To compute the Hilbert-Schmidt degree of non-Gaussianity we 
have to use in Eq. ( 12. lip the overlap ( 14. 7ft . the purity ( 14. 9p of the damped state, and 
the purity of the reference Gaussian state, 



Tr 



P ( a M) 



(t) 



(4.10) 



2{N(t)) + 1 

In Fig.3 we show the time evolution of the Hilbert-Schmidt degree of non-Gaussianity for 
the same sets of states as in Figs.l and 2. The agreement between the Hilbert-Schmidt 
degree of non-Gaussianity and the other two distance-type measures we here discuss is 
more accurate during the interaction with noisier reservoirs, as shown in Figs. 2 and 3. 



4-3. Gaussification and entropy production 

In the following we shall take advantage of having evaluated the von Neumann entropy 
for both the time- dependent state p^ M \t) and the reference Gaussian one Pa{t). In 
the latter case, the field-reservoir interaction preserves the feature of p^it) of being 
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yt yt 

Figure 3. Time evolution of the Hilbert-Schmidt degree of non-Gaussianity for contact 
with a heat bath when Tlr — 0.1 (left plot) and ur — 5 (right plot). The input 
states are A/-photon-added thermal states with n = 1 (left) and n = 0.5 (right), and 
M = 1, 3, 5, 10 (from bottom to top) 

a thermal state. Its mean photon number (N(t)), Eq. (14. 5p . varies monotonically from 
the input value (N(0)) = M(n + 1) + n to the reservoir mean occupancy n^. Therefore, 
the von Neumann entropy S(pq (t)) is a monotonic function of time between the 
corresponding imposed limits at t = and t — > oo. Good examples of this monotonic 
evolution are shown in the left-hand side of Fig. 4 for a noisy reservoir with = 5. The 




yt yt 

Figure 4. Entropy production during the gaussification presented in Fig. 2, namely: 
fin = 5, n = 0.5, and M = 1, 3, 5, 10 (from bottom to top) 

input mean occupancy (iV(0)) for the states with n = 0.5 and various M are placed in 
Table 1 together with the corresponding values of the entropy for the reference Gaussian 
states Pq ^ and the given ones . To calculate the latter ones we have used the density 
matrix, Eq. (14. 4p . Notice that in all the chosen cases the input entropy of 

p{M) is 

smaller 

than that of the reservoir. The von Neumann entropy of the damped state is plotted in 
the right-hand side of Fig.4 for the same parameters. 

We point out the existence of a short-time maximum of entropy when (iV(0)) > hr 
(cases M — 5, 10 as seen in Table 1). A similar behaviour was previously reported 
for some pure input states: a number state in Ref.[36] and an even coherent state in 
Ref.[37J. It was found analitically that a maximum in the evolution of the 2-entropy of 
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M 


(N(0)) 






1 


2 


1.909 


1.372 


3 


5 


2.703 


1.824 


5 


8 


3.140 


2.078 


10 


15.5 


3.770 


2.429 



Table 1. Entropic characterization of the input states whose gaussification is analyzed 
in Figs. 2 and 3, namely having n = 0.5. The thermal mean occupancy of the reservoir 
is ur = 5. 

the attenuated mode of these pure non-classical state exists provided that their initial 
mean photon number exceeds the mean occupancy of the reservoir. 

5. Discussion and conclusions 

All plots showing the time development of 5p and 5 re in Figs.l and 2 have a similar 
aspect which indicates a good consistency of these two measures coming from their 
monotonicity under the damping map. We expect they induce the same ordering of 
non-Gaussianity. This property is not shared by the Hilbert-Schmidt degree of non- 
Gaussianity and we have to remark the different concavity of the left plots in Fig.3. For 
weak damping this measure could not induce the same ordering of the non-Gaussianity 
as the relative entropy or Bures distance. In Fig.5 (left) we have plotted the entropic 
measure as a function of the fidelity-based degree for the set of damped states with 
ur = 0.1, n = 1 and different values of M. The two measures show here a good 
consistency by preserving the same ordering of non-Gaussianity. In Fig.5 also plotted 
are 5 re an d 6p as functions of 5rs f° r the same sets of states as in the left-hand side 
panel. We notice that both of them have the same concavity but they lack consistency 
for some values of the parameters, especially at small times. 




<5f £>hs <5hs 



Figure 5. Displaying a full consistency between 5re and 5p (left panel) and a partial 
one between these two measures and the Hilbert-Schmidt one for input M-photon 
added thermal states (n = 1) evolving in contact with a thermal bath when hr =0.1. 
We used M = 1,3, 5, 10 (from bottom to top). 

To conclude, in this work we have investigated the evolution of non-Gaussian states 
during the interaction of the field mode with a thermal reservoir resulting in two 
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processes: loss of non-Gaussianity shown by some recently introduced distance-type 
measures, and a non-monotonic time development of the von Neumann entropy. For 
the specific states we were interested in, we have found that decaying of non-Gaussianity 
was consistently described by the three distance-type measures used in this paper. Even 
with some alterations of ordering, the fidelity-based degree, the Hilbert-Schmidt one and 
the entropic measure evolve monotonically, as expected for any good measure of non- 
Gaussianity [19]. On the contrary, the behaviour of the von Neumann entropy of the 
damped state depends on the relation between the input mean occupancy (iV(0)) of 
the state and the thermal mean occupancy % of the reservoir. States with input mean 
occupancies greater than ur, but with input entropy not exceeding that of the reservoir 
are found to present a transient mixing enhancement, visible as a maximum in the 
evolution of their entropy. 
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Appendix A. Two power series involving Gauss hypergeometric functions 

A Gauss hypergeometric function is the sum of the corresponding hypergeometric series, 
2 FM^c,z):=^^f^- r (|*| <1), (A.l) 

where (a) n := T(a + n)/T(a) is Pochhammer's symbol. This definition is extended by 
analytic continuation. The Legendre polynomial of degree M can be expressed in terms 
of a Gauss hypergeometric function: 

Pm{z) = 2*1 —M] 1; Z -^j , (M = 0, 1, 2, 3, ...).(A.2) 

Let us consider a well-known generating function of a class of hypergeometric 
polynomials [37] : 

oo 

G {u,z) :=^u l 2 F 1 (-l,a;l;z) = {l-u) a -\l-u+uz)~ a , \u\ < 1.(A.3) 

1=0 

In the body of the paper we need to evaluate the sum 



G n (u, z):=f:( l n )u l 2 F X {-1, a; !;*) = ! „« z) 

l=n \ / 



(A.4) 



We have performed the sum (1A.4|) in two ways. First, we have applied an elegant method 
described in the recent paper [48J to find the rath-order derivative of such a generating 
function using Cauchy's integral formula. Second, because the r.h.s. of Eq. ( 1A.3I) has a 
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simple structure, the nth-order derivative can be put in a closed form after some routine 
algebra. With both methods we get: 

GJu,z) =u n (l-u) a ~ n - 1 (l-u + uz)- a 2 F 1 ( -n,a;l; j .(A.5) 

\ 1 — u + uz J 

The known summation of a power series whose coefficients are proportional to the 
product of two Gauss hypergeometric polynomials [UJ is also employed in Sec. 4: 



ra=0 



ST 2 F 1 (-n, 6; c; ^(-n, 6'; c; z') 



+ 
x 2 Fi 



-c+6+b' 



(1+f -fz)" b (l+f -f^') 



6, 6'; c; 



— fzz 



(1 + u — vz){l + v — vz') 
In our particular case b' — b, c — 1, z' — z, so that Eq. (1A.6|) becomes 



(A.6) 



J2 v n [ 2 Fi(-n, 6; 1; ^)] 2 = (1 - u)- 1+26 (l - « + vz)' 



2b 



n=0 



X 2^1 



6,6;1; 



t>2 



'1 — f + t> z) 



(A.7) 
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